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Abstract 

We calculate the lowest order hadronic cross sections for producing colored heavy 
quark-antiquark pairs in L = 5" = and L = S = 1 configurations. Such QQ[^Sq^^] and 
QQl'^P^f''^] states hadronize into i/jq quarkonia at the same order in the NRQCD velocity 
expansion as previously considered (5(5['^5'^^''] pairs. Their contributions to prompt Psi and 
Upsilon production at the Tevatron bring the shapes of theoretical transverse momentum 
distributions into line with recent CDF measurements. We find that the best fit values for 
the linear combinations of QQ[^Sq^^] and QQ[^Pj^^] long distance matrix elements which 
can be extracted from the data are generally consistent with NRQCD scaling rules. 



4/96 

^ Work supported in part by a DuBridge Fellowship and by the U.S. Dept. of Energy under 
DOE Grant no. DE-FG03-92-ER40701. 

^ Work supported in part by the U.S. Dept. of Energy under DOE Grant no. DE-FG03-92- 
ER40701. 



1. Introduction 

The study of quarkonia has yielded valuable insight into the nature of the strong inter- 
action ever since the discovery of the J /ij) resonance in 1974. During the past two decades, 
QQ bound states have provided useful laboratories for probing both perturbative and 
nonperturbative aspects of QCD. Recently, investigations of charmonia and bottomonia 
systems have uncovered some striking surprises. Orders of magnitude disagreements have 
been found between old predictions and new measurements of Psi and Upsilon production 
at several collider facilities. These large disparities have called into question the simplest 
model descriptions of quarkonia and stimulated the development of a new paradigm for 
treating heavy quark-antiquark systems based upon QCD. Although much theoretical and 
experimental work remains to be done before a truly consistent picture of quarkonia is 
established, ongoing studies of these heavy mesons are leading to a better understanding 
of some basic aspects of strong interaction physics. 

Quarkonia bound states are qualitatively different from most other hadrons since 
they are inherently nonrelativistic. The physics of quarkonia consequently involves several 
energy scales which are separated by the small velocity v of the heavy constituents inside 
QQ bound states. The most important scales are set by the mass Mq, momentum MqV 
and kinetic energy Mqv'^ of the heavy quark and antiquark. In order to keep track of this 
scale hierarchy, an effective field theory called Nonrelativistic Quantum Chromodynamics 
(NRQCD) has been established [0. This effective theory for QQ bound states shares several 
similarities with the Heavy Quark Effective Theory (HQET) which describes the low energy 
QCD structure of heavy-light Qq mesons. For example, NRQCD is based upon a double 
power series expansion in the strong interaction fine structure constant ctg = gl/^'n: and 
the velocity parameter v ~ 1/ log Mq which is similar to the HQET's double expansion 
in Us and 1/Mq. Both theories also incorporate approximate spin symmetry relations 
which constrain various multiplet structures and transition rates. But most importantly, 
NRQCD systematizes one's understanding of charmonia and bottomonia just as the HQET 
methodically organizes the physics of D and B mesons. 

Quarkonia are described within the NRQCD framework in terms of Fock state decom- 
positions. The wavefunction of an S-wave orthoquarkonium vector meson schematically 
looks like 

I ^q) = 0(1) I QQ p^f ) + 0{v) I QQ ['Pf] g) 

+ 0{v') I QQ g) + 0{v') I QQ p^^'^] gg) + 0{v') \ QQ fD^j^'^] gg) + ■ ■ ■ . 

(1.1) 
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The spin, orbital and total angular momentum quantum numbers of the QQ pairs in each 
Fock component are indicated within the square brackets in spectroscopic notation, while 
the pairs' color assignments are specified by singlet or octet superscripts. The order in 
the velocity expansion at which each of these Fock states participates in i/jq annihilation 
or creation processes is governed by simple NRQCD counting rules [^j. For instance, 
suppose a heavy quark and antiquark are produced almost on-shell with nearly parallel 
3-momenta in some high energy reaction. The low energy hadronization of this pair into 
a physical i/^q bound state takes place at 0{v^) if it has the same angular momentum and 
color quantum numbers as those displayed in the first Fock component of (|1 . 1| ) . The long 
distance evolution of all other QQ pairs generated at short distance scales into i/jq mesons 
occurs at higher orders in the velocity expansion. 

If the relative importance of various quarkonia production channels depended solely 
upon the order in v at which pairs hadronize into physical bound states, those modes 
which proceed through the leading Fock components in quarkonia wavefunction decompo- 
sitions would generally be dominant. This assumption coincides with the basic tenet of 
the so-called color-singlet model Quarkonium production is presumed in this model 

to be mediated by parton reactions that generate colorless heavy quark-antiquark pairs 
with the same quantum numbers as the mesons into which they nonperturbatively evolve. 
Transverse momentum distributions calculated within this picture badly underestimate 
experimental observations for p± > 2Mq. The breakdown of the color-singlet model stems 
from its neglect of all high energy processes that create QQ pairs with quantum numbers 
different from those of the final state meson. In particular, it overlooks short distance 
contributions to quarkonia cross sections from intermediate color-octet states which can 
be orders of magnitude larger than those from color-singlet pairs. Even if the long distance 
hadronization of the former is suppressed by several powers of v compared to the latter, 
the color-octet components of quarkonia distributions can dominate overall. 

In a previous paper 0], we examined the contributions to hadron collider quarkonia 
cross sections from processes which create colored heavy quark-antiquark pairs in L = 0, 
S = 1 configurations at short distance scales. Such QQ['^s[^^] pairs may emit or absorb 
a soft chromoelectric dipole gluon and evolve at 0{v^) into Xqj mesons. These P-wave 
quarkonia can later radiatively decay to lower S-wave V'q levels. Alternatively, QQ[^s[^^] 
pairs may undergo double chromoelectric dipole transitions and directly hadronize at 0(f ^) 
into i/jq mesons. The first channel formally represents the dominant color-octet production 
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mechanism. But since Xqj ~^ V'q + 7 branching ratios in both the charmonia and bot- 
tomonia sectors are numerically comparable to f ^ ~ 0.23 and f ^ ~ 0.08, the second mode 
is phenomenologically important as well. Consistency then requires that 0(f ^) color-octet 
contributions to i/jq production from colored L = S = and L = S = 1 states also 
be considered. Therefore, we calculate in this article the i/'q cross sections arising from 
QQ[^Sq^^] and QQ[^Pj^^] intermediate channels. As we shall see, inclusion of a fully con- 
sistent set of color-octet quarkonia cross sections brings theoretical predictions into line 
with experimental observations of prompt Psi and Upsilon production at Fermilab. 

Our paper is organized as follows. In section 2, we compute the amplitudes for gg — > 
QQ[^Sq^^] and gg — > QQ[^P^f'^] scattering which mediate if^Q production at small transverse 
momenta. We then utilize these amplitudes in section 3 to determine the QQ[^Sq^^] and 
QQ[^Pj^^] contributions to i/'q differential cross sections in 2 — > 2 collisions. In section 4, 
we fit color-octet quarkonia distributions to recent Psi and Upsilon measurements and 
determine improved numerical values for several NRQCD matrix elements. Finally, we 
summarize our findings in section 5 and close with some thoughts on the implications of 
our results for quarkonia production in several different experimental settings. 

2. tpQ production in 2 — » 1 collisions 

Color-octet quarkonia production starts at 0{a'^) with the scattering processes 13'+^ 
Q + Q and g + g ^ Q + Q- The Feynman diagrams which mediate these reactions are 
illustrated in fig. |I| and fig. ^ In the first quark scattering channel, the heavy quark- 
antiquark pair appearing in the final state must share the same quantum numbers as its 
intermediate virtual gluon progenitor. Angular momentum, parity and charge conjugation 
conservation restrict the spin, orbital and total angular momentum quantum numbers of 
the QQ pair toL = 0or2, 5' = 1 and J = 1. In the second gluon scattering process, the 
allowed values for L, S and J are not so tightly prescribed. In ref. |]^, we examined the 
important special case where gluon fusion produces QQ[^s[^^] pairs. We now generalize 
our earlier analysis and consider the 0{a'j.) formation of other pairs with different sets of 
quantum numbers. 
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To begin, we write down the on-shell scattering amplitude 



X U 



8M2 



1 p 

fabciTcTj [g^iuih - i>i) + 2(pi^7ax - P2m7!^)] |^(-^ - 9; S2) 

(2.1) 

that corresponds to the sum of the three graphs displayed in fig. ||. The reduced amplitude 
which describes the creation of a QQ pair in a particular angular momentum and color 
configuration is obtained from this expression by applying a series of projection operations: 



Lz,Sz si,S2 i,j 



p p 

X A ( ga{pi)9b{p2) ^Q'i^ + q; si)Qj{— - q\ S2, 



(2.2) 

Several points about this projection formula should be noted. Firstly, a QQ pair has 
negligible overlap with a nonrelativistic quarkonium bound state unless the relative mo- 
mentum q between the heavy quark and antiquark is small compared to their combined 
momentum P . We have therefore incorporated a delta function into (|2.2|) which restricts 
the triple integral over q to the two-dimensional surface defined by \q \ = a/ Mq^ where 
qO Mv'^ <^ M = 2Mq. The squared invariant mass of the QQ pair thus equals = 
up to small relativistic corrections. Inclusion of the delta function also properly converts 
the mass dimension of the 2-^2 scattering amplitude into that for a 2 1 reaction. 
Secondly, the angular integration over the spherical harmonic projects out the pair's spec- 
ified partial wave. The sums over the Clebsch-Gordan coefficients similarly project out 
the spin and total angular momentum of the QQ object. Finally, the sum over the SU{3) 
coefficients 

(3i;3j|l) = 5?/v^ 

(2.3) 

{3i;3j\8c) = V2{T,)i 

combines together the color quantum numbers of the quark and antiquark into either a 
singlet or octet configuration. 



Inserting the gluon scattering amplitude into the projection formula in ( |2.2|) , 

we can readily calculate the reduced amplitude for two gluons to fuse into an arbitrary 
color-octet combination. We list below the formation amplitudes for QQ[^s[^^], QQ[^Sq^^] 
and QQ[^Pj^^] pairs which all hadronize into i/jq bound states at O(f^) in the velocity 
expansion: 



^ 2/2 [(^^2 - Pi)agf„y + 2{pi„gf,^ - P2^^gua)] (2.4a) 

pi + pi — -' 

M i \ i \ nniiQi^)! \ ^ /47rM ^ Me^r^uPiPle>'{pi)e'' {P2) , . 
Aiga{pi)gb{P2) ^ QQ[ S'^ %) = -^^^J—^gj^,^ p2^p2_M2 (2-4^) 

A{ga{Pi)gb{P2)^QWPj%) = E {^L,;lS,\JJ,)e'^{L,re^{S,r 

Lz,Sz 

Me''{pi)e''{p2) 

X 



2 , 2 7Tr\9a^iPl,,iP2 - Pl)p + gcuP2^i{Pl - P2) (3 

Pi + P2 -'^^ ^ 



+ {pI - Pi - M'^)gcvgp^, + [pI -pI- M'^)g^^gp^ 

+ A^2 ^2 2 [d^i^PlpiPl - P2)c + g^luP20{p2 - Pl)a] 

M"^ — Pi — P2 



^ m^+pI+pI 



_p2_ p2 [9t^pPlu{P2 - Pl)a + gu/3P2^,{Pl - ^2)0] }■ (2.4c) 

As required by gauge invariance, these expressions vanish when Pi=P2 = and e^{pi) —>■ 
Pi or e^{j>2) P2- The general P-wave result listed in ( |2.4c| ) may be further reduced by 
employing the Clebsch-Gordan identities 

Lz ,Sz 

J2 {^L,- lSz\lJz)e-{L,Ye^{Szr = -^'"^''^1^^^^-)* (2-5) 

Lz ,Sz 

J2 {lLz;lSz\2Jz)e"{Lzre^{Szr = e'^^iJzr. 

Lz,Sz 

We then find that gg QQl^P^^^] as well as gg QQ[^s[^^] scattering vanishes when 
both incident gluons go on-shell 

The projection formula in (^T^) can obviously be generalized to other parton channels 
besides gg ^ QQi^^+^L^j'^^]. We may insert any QCD amplitude which has a heavy 



quark and antiquark appearing in the final state and project out a reduced color-singlet or 
color-octet expression. Applying this general technique to the qq QQ scattering process 
pictured in fig. |I|, we find 



A{q{pi)q{p2) - QQ?S[\) = ^^^^9XP2haTau{p,)e-{S^r. (2.6) 

In ref. , we performed a less general projection operation which did not explicitly involve 
a relative momentum integration nor a partial wave decomposition of A{qq —>■ QQ). The 
corresponding reduced amplitude listed in eqn. (2.10) of ref. [§] thus differs from the result 
displayed above by an overall multiplicative factor. 

The squares of 2 — > 1 amplitudes enter into the differential cross section for heavy 
pair production 

da (a{p^)b{p2) - QQi'^^+'L^j^'^m) = ^ ^ U(a6 ^ QQ^^+^lJ/'^)]) ' d^^ip^+p2; P) 
\ / ^PrP2 — I 

(2.7) 

where the barred summation symbol indicates that initial (final) spins and colors are 
averaged (summed) and d^i denotes a one-body phase space factor. High and low energy 
effects are intertwined in this expression. In order to disentangle it, we follow ref. |^ and 
match the integrated cross section onto the product of a short distance coefficient and a 
long distance NRQCD matrix element: 

a (ab - QQ[''^'L^J''h) = §S x {0\O?§C'+'Lj)\0). (2.8) 

Mq 

The general structure of the operator whose vacuum-to-vacuum matrix element appears 
on the RHS of this matching condition looks like 

Of^C'^'Lj) = x'Ki: {Y\QW^'Lf'^]) {qW'^'L^}^'^]\\ ^'^x (2.9) 



where and x represent two-component Pauli spinor fields and the matrix K denotes a 
product of color, spin and covariant derivative factors. The intermediate quark-antiquark 
state sandwiched in the middle 

73„ 7f2. 



Lz,Sz si,S2 i,j 

xYlL^{q){^si;^S2\SSz){LLz;SSz\JJz){3i;3j\hS) Q\q; si)Q^{-q; S2) 



(2.10) 
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is defined in the NRQCD effective theory in the same way as in fuU QCD. As a result, 
any arbitrariness in the definitions of the heavy pair production cross section and NRQCD 
matrix element cancels out of their ratio. The short distance coefficient appearing on the 
RHS of ( p78|) is convention independent. 

All information related to the hard scattering process which creates the QQ pair is en- 
coded within Cshort- This same coefficient enters into the physical quarkonium production 
cross section 

a(a6-^QQp^+^4^'«)]-^,+x) =^^x (Oiotir^^^ (2.11) 

On the other hand, the accompanying long distance matrix element which specifies the 
probability that a QQ['^^~^^L^j'^^] pair hadronizes into a i/jq bound state is completely 
different from its counterpart in ( p.8|) . The operator 

Ofj C'^'Lj) = x^Kij E I^'^ + ^) ^^'^ + ^1 ) ^^^^ (2.12) 

\mj X J 

has one unit greater mass dimension than 0^g{'^^~^^Lj) as can be verified by comparing 
the dimensions of heavy intermediate pair and nonrelativistically normalized i/'q states. 
The inverse powers of Mq in cross section equations (|2.8|) and (|2.11|) consequently differ 



by unity. The nonperturbative matrix element (O|0f § (^'^+^Lj)|0) also cannot readily be 
calculated within NRQCD unlike its perturbative {0\Oi^ (^^~^^ Lj)\0) counterpart. Simple 
multiplicity relations such as 

(010(^8(3^010) = 3 (0|Of8(i5o)|0) ^ ^ 

2.13 

mO^^si'Pjm = (2J + l)(0|Of8(=^Po)|0) 

are obeyed exactly by the latter and approximately by the former. But the color-factor 
relation 

]\t2 _ 1 

{0\OiC'+'Lj)\0) = -^^{0\O^C'+'Lj)\0) (2.14) 

which holds for H = QQ certainly does not apply when H = i/^q. Numerical values for 
{0\Of'g(^^~^^Lj)\0) matrix elements must be extracted either from experimental data or 
lattice calculations. 

In order to clarify the meaning of these NRQCD matching ideas, we explicitly evaluate 
the matching conditions specified in eqns. ( p.8|) and (|2.11| ) for one simple example. We 



consider the gluon fusion formation of an ryg pseudoscalar meson through an intermediate 
QQI^Sq'^^] pair. A straightforward computation yields the color-singlet cross section 



384712 qOg 
and matrix element 

{o\o?^CSom = (o QQl'si,'^]) (qqP^J') 



(2.15) 



(2.16) 



We derived this last result within the NRQCD effective theory by decomposing the Pauli 
fields 



X 



(2.17) 



^ J (27r 

in terms of two-component spinors normalized according to 

2 2 

X^^a(p, s)4(p, s) = ^r]c,{p, s) r]l{p, s) = 6ap 



(2.18) 



s=l 



s = l 



and single fermion creation and annihilation operators which satisfy the nonrelativistic 
anticommutation relations 



{b\p;s),b]{p';s')} = {c\p;s),cUp';s')} = {27rfd;d,,,d^'Hp-p'). 



(2.19) 



Taking the ratio of ( |2.15| ) and ( p.l6| ) , we deduce the short distance coefficient in matching 
condition 



C{gg - QQ[^4'^])short _ 1 Tr^a^ ^ 



M, 



d-4 

Q 



9 M^s 



(2.20) 



and the gluon fusion cross section in matching condition ( |2.11|) : 

a (^7^7 - QQi'si'^] - r?^) = (i - ^) {o\oresom. (2.21) 



If we recall the relation between the NRQCD matrix element and squared rjg wavefunction 
at the origin [^j 

(2.22) 



{0\OrCSom = ^R{0)\l + Oiv'))., 



we see that our result is consistent with the 0{a'^) cross section 



a (^gg QQi^Sl,'^] ^ 



) 




(2.23) 



3 M^s 



previously reported in the literature 0. 

Working in a similar fashion, we can decompose any color-singlet or color-octet cross 
section into products of short and long distance factors. We tabulate in the Appendix all 
0(Q!g) short distance squared amplitudes for 2 ^ 1 color-octet reactions which yield i/jq 
bound states at 0{v''') in the NRQCD velocity expansion. The corresponding long distance 
factors are simply given by appropriate NRQCD matrix elements for specific production 
channels. For example, the total squared amplitude for gg — > QQ[^Sq^'^] ^ipQ scattering 
equals the product of the process-independent high energy expression listed in eqn. (A. 2a) 
and the process-specific low energy matrix element (0|Og '^(^5'o)|0): 



Color-octet pair production in 2 — 1 collisions could represent an important source 
of quarkonia in fixed target experiments, and its impact needs to be studied. But before 
definite predictions can be made, numerical values for color-octet matrix elements must 
be known. Therefore, we now turn to consider quarkonia production at hadron colliders 
where we can use experimental data to determine these matrix element values. 

3. tpQ production in 2 — » 2 collisions 

In order to be experimentally detectable, quarkonia must be created at collider facil- 
ities with nonvanishing transverse momenta so that they are not lost down the beampipe. 
Hadrons resulting from 2 — > 1 scattering processes typically have small p±_ comparable 
to the QCD scale. The production of quarkonia with nonnegligible transverse momenta 
therefore mainly proceeds through 2^2 collisions. Such reactions start at 0(q:^) via the 



parton channels qq QQ[^^+^L^j'\, gq QQ[^^+^L^j'% and gg QQ[^^+^L^]-'^^]g. 



In ref. [||, we calculated the differential production cross sections for L = 0, = 1 color- 
octet pairs. In this section, we generalize our previous results and consider the formation 
of colored L = S = and L = S = 1 pairs. 

The Feynman diagrams which mediate quarkonia production in these color-octet chan- 
nels are illustrated in fig. |^. The shaded circles appearing in the figure represent the 



J2\A{gg^QQ['si,'^]^i;,,) 



5(47r«,)^ 
192M 



{oiof'Csom. 



(2.24) 
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gg -> QQI^S'-Q^ and gg QQ[^Pf^] amplitudes in eqns. ( CT) and (|^. The qq i^gg 
and gq i/^Qq diagrams pictured in figs. 3a and 3b can readily be squared using standard 
spinor summation techniques. On the other hand, conventional evaluation of the gluon 
channel graphs in fig. 3c represents a formidable computational task. It is therefore ad- 
visable to find a more tractable method for calculating the color-octet contributions to 
99 — ^ i^Q9 scattering. 

We adopt a simple helicity amplitude technique to sum and square the gluon graphs 
in fig. 3c. We first choose the following explicit representations for the gluon momenta and 
polarization vectors shown in the figure: 



Pi 

P2 
P4 



(1,0,0,1) 

(1,0,0,-1) 



1 



V2 



(0,1, z,0) 



V2 



(0,1, -i,0) 



(3.1) 



(l, 0, sin6', — cos 6*) 



V2 



(0, ±1, i cos 6*, i sin6'). 



We next boost the the heavy pair's four-momentum from its primed rest frame to the 
unprimed lab frame: 



p^ = (M,0,0,0)^P3 



-ITT'"' 



sin 6*, 



cos 9 



2V§ ' 2V§ 

We also Lorentz transform the rest frame polarization vectors and tensors 



1 

—i 
\ / 



^0 ON^ 

-i 

-i 

Vo i/ 



(3.2) 



[e 



I (/l = 0)^ 






\\) 



/ {h=±\)Y 




(3.3) 



(4 



1 

i 

Vo/ 













°\ 


1 





1 







2 





T« 


-1 







Vo 








0/ 



of J = 1 and J = 2 QQ pairs. Given these explicit representations, it is easy to work 
out all possible scalar contractions and express the answers in terms of the Mandelstam 
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invariants s, t = -{s - M'^){1 - cos6')/2 and u = -{§ - M'^){1 + cos6')/2. The gluon 
channel amphtudes are functions of these Lorentz invariant dot products. 

Using the high energy physics package FEYNC ALC , we calculated each individual 
helicity amplitude for gg QQ[^SQ^^]g and gg QQ[^Pj^^]9 scattering. Parity and 
crossing symmetry relations between different helicity amplitudes provided valuable checks 
on our Mathematica code. Since separate helicity amplitudes do not interfere, the total 
squared amplitude simply equals the sum of the squared helicity amplitudes. The final 
results for the gg QQ[^^~^^L^p]g channel are displayed in the Appendix alongside those 
for the qq — > QQ[^^~^^L[f '^]g and gq — QQ['^^~^^L[^'^]q modes. For completeness, we include 
in this list the QQ['^S[^^] squared amplitudes which we calculated in ref. [If. 

The products of short distance color-octet squared amplitudes and long distance 
NRQCD matrix elements enter into the partonic cross section 



(3.4) 

After folding in distribution functions fa/Ai^a) and fb/Bi^b) that specify the probabilities 
of finding partons a and b inside hadrons A and B carrying momentum fractions Xa and 
Xb, we obtain the hadronic cross section 



dy3dy4dp± ' ' dt 

abc 

(3.5) 

which is a function of the tpQ and recoiling jet rapidities and 1/4 and their common 
transverse momentum p^. With this hadronic distribution in hand, we can determine 
color-octet contributions to -i/^g production in any hadronic process. We apply it to the 
study of charmonia and bottomonia at Fermilab in the following section. 



4. Psi and Upsilon production at the Tevatron 

During the past few years, striking disparities have arisen between old predictions and 
new measurements of J/V', V'' and T production at the Tevatron. The CDF collaboration 
has detected these heavy mesons at rates which exceed theoretical expectations based upon 
the color-singlet model by orders of magnitude [1T0|-[T2[| . In ref. , we examined the impact 
of cc['^S'j^'*] and intermediate states upon Psi and Upsilon production. Since 
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numerical values for most NRQCD color-octet matrix elements were unknown, we simply 
fitted the magnitudes of da/dp±{pp ^ i/jq + X) octet cross sections to the CDF data. We 
found that including the QQf'^S'j;^'*] channel significantly diminished discrepancies between 
the shapes of the predicted and measured transverse momentum distributions. We now 
update our earlier analysis and incorporate cross section contributions from QQ[^Sq^^] and 
QQ[^Pj^^] pairs. As we shall see, the fully consistent O(f^) set of color-octet differential 
cross sections yields substantially improved fits to the data. 
We first plot in fig. || the ratio 

Rip±) = (4.1) 

where we temporarily set (0^'^(3Po)) = M^iOt'^CSo)) for comparison purposes. The 
solid curve's nearly constant value R{p±) ~ 3 for pj^ > 5 GeV indicates that the shapes 
of the cc[^Sq^'^] and cc[^Pj^'^] differential cross sections are practically identical in the 
charmonia sector. As a result, all fits for the NRQCD matrix elements in these color- 
octet channels become degenerate when performed over the transverse momentum range 
5 GeV < p± < 20 GeV where J/V' and differential cross sections have been measured. 
We consequently can only extract the linear combination {Ot{^Po))/M^ + {0^{^So))/S 
along with {Og{^Si)) from the CDF data. In the bottomonia sector, the shapes of the 
bb[^SQ^^] and bb[^Pj^^] distributions are not exactly the same throughout the < p± < 
15 GeV interval where Upsilon data exists. As indicated by the dot-dashed curve in 
fig. ^ R{p±) varies around 5 over this transverse momentum range. Yet the differences 
in shape between the bb[^Pj^^] and 66[^S'q^'*] contributions to the total Upsilon differential 
cross section are not sufficiently great so that a full three-parameter color-octet matrix 
element fit can be reliably performed. So we will simply determine estimates for the linear 
combination {Oj{^Po))/M^ + {Oj{^So))/5 along with {Oj{^Si)). 

Our new fits to prompt charmonia production at the Tevatron within the pseudora- 
pidity interval |r7| < 0.6 are illustrated in fig. || and fig. ^ All contributions from B meson 
decay have been removed from the data sets displayed in these figures, and radiative Xcj 



^ The diflFerential cross sections which enter into results displayed in fig. ^ and ah subsequent 
figures were calculated using the MRSDO parton distribution functions evaluated at the renor- 
malization scale fi = a/Pj^^K/W^. 



12 



decay feeddown to the J/?/' differential cross section has been separated out as well. The 
dashed curves depict the direct color-singlet production predictions based upon the charm 
quark mass value Mc = 1.48 GeV and the Buchmiiller-Tye charmonium wave functions 



at the origin tabulated in ref. [|T3[. The dot-dashed and dotted curves illustrate the best 
fits for the ccpS'^^''] and combined cc[^Pf''^] plus cc[^Sq^^] channels. The solid curves show 
the sums of the color-singlet and color-octet components and represent the total predicted 
differential cross sections. 

Following the interpolation procedure described in ref. we have included leading 
log corrections into the ccf^S"!^^^] differential cross sections so that they approach Altarelli- 
Parisi improved gluon fragmentation distributions for pj_ ^ Mc- In the large transverse 
momentum limit, gluon fragmentation represents the dominant source of prompt charmo- 
nia ||T3|-pil[1. This asymptotic behavior can be seen in the dotdashed cc['^5'|^^''] curves of 
figs. 5 and 6. But throughout the < < 20 GeV region, they are not overwhelm- 
ingly larger than the combined cc["^pj^'*] and cc[^Sq^'^] components whose contributions to 
prompt charmonia production are sizable. Inclusion of the latter color-octet channels into 
the total differential cross sections yields theoretical t/j' and J/ tj} distributions which fit the 
data quite well. Their respective x^/NDOF = 0.5 and x^/NDOF = 0.9 figures-of-merit 
are nice and small. 

In fig. ^ we plot the transverse momentum distribution of J/ij) mesons which result 
from radiative Xcj decay. The dashed curve in the figure shows the color-singlet Xcj 
differential cross section multiplied by Br(xcj J /i^ + l) and summed over J = 0, 1 
and 2. The dot-dashed curve illustrates the cc['^<S'^^''] channel contribution. The solid 
curve corresponds to their sum and represents the total 0{v^) cross section prediction. 
As indicated by its poor x^/NDOF = 2.3 value, this solid line does not fit the data well. 
We believe that a better match could be achieved if subleading color-octet contributions 
were included. The first subdominant corrections enter at 0(f ^) in the NRQCD velocity 
expansion from the long distance evolution of QQ'^pf'\, QQpDj^'*] and QQ\^P^^] pairs 
into Xqj bound states. Since short distance production cross sections for the latter two 
pairs have not yet been calculated, we cannot legitimately include into fig. |^ subleading 
contributions from the first pair which we have computed. 

We turn now to the bottomonium sector and consider Upsilon production at the Teva- 
tron within the rapidity interval \y\ < 0.4. Our new fits to CDF T(15') and T{2S) data 
are displayed in fig. || and fig. |^. No separation between prompt and delayed Upsilon 
sources has been experimentally performed. The dashed curves in the figures therefore 
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include both direct T production and radiative feeddown from Xhj states. These color- 
singlet distributions are based upon the bottom quark mass value Mh = 4.88 GeV and 
the Buchmiiller-Tye bottomonia wavefunctions at the origin tabulated in ref [0. The 
dot-dashed and dotted curves illustrate the hh[^sf\ and combined hh[^Pf\ plus hh[^s'i\ 
fits. The solid curves equal the sums of the color-singlet and color-octet contributions and 
represent the total T differential cross sections. As we previously discussed in ref. [§] , the 
color-singlet and color-octet distributions are corrupted at very small transverse momenta 
by collinear divergences which should be factored into incident parton distribution func- 
tions. Soft gluon effects also need to be resummed before the cross section turnover which 
is evident in fig. H can be properly described. Since we have not incorporated these effects, 
our cross section predictions are not trustworthy at low p_^. We therefore exclude points 
in fig. H and fig. ^ with p±_ < 3.5 GeV from our fits. We then find x^/NDOF = 0.3 and 
X^/NDOF = 0.9 for the remaining points in these figures. 

NRQCD power counting rules provide useful consistency checks on all our fits. We list 
in tables I and II the numerical values for color-octet matrix elements which we extracted 



Color-Octet 
Matrix Element 


Numerical 
Value (GeV^) 


NRQCD 
Scaling Order 


(0108^/^(3^1)10) 


(6.6 ± 2.1) X 10-2 


MM 




(9.8 ± 1.3) X 10-3 


MM 




(4.6 ± 1.0) X 10-3 


MM 




(5.9 ± 1.9) X 10-3 


MM 


^0|Og^"(i^)(3^,)|0) 


(4.2 ± 1.3) X 10-^ 


MM 




(4.1 ±0.9) X 10-3 


MM 


(0|Of^('^)(3^i)|0) 


(3.2 ± 1.9) X 10-1 


MM 



Table I. Color-octet matrix elements 
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Color- Octet Matrix Element 
Linear Combination 


Numerical 
Value (GeV^) 


NRQCD 
Scaling Order 


{0\Oi/^{^Pom (0108^/^1^0)10) 
M2 3 


(2.2 ±0.5) X 10-2 




mofePom {o\ofcsom 

3 


(5.9 ± 1.9) X 10-3 


MM 


5 


(7.9 ± 10.0) X 10-3 


MM 


5 


(9.1 ± 7.2) X 10-3 


MM 



Table II. Color-octet matrix element linear combinations 



from the data along with their scaling dependence upon the heavy quark mass Mq and 
velocity Vq. The values for all the charmonia matrix elements were derived directly from the 
CDF J /ij) and tp' data. On the other hand, insufficient experimental information exists to 
independently extract (0|Oj^''^^(35i)|0) and (0|O^"^''-^^(3^i)|0) in the bottomonia sector. 
We therefore determined the latter from the Upsilon data after having scaled up the former 
from the corresponding Psi color-octet matrix elements using NRQCD power counting 



rules. The remaining (0|O 



T(nS).3 



Po)|0)/M2 + (0|0, 



r{uS)n 



5'o)|0)/5 linear combinations 



were obtained directly from the bottomonia cross section data. 

The error bars listed in tables I and II are statistical and do not reflect systematic 
uncertainties in heavy quark masses, color-singlet radial wavefunctions, parton distribu- 
tion functions and next-to-leading order corrections. The magnitudes of all these dif- 
ferent sources of uncertainty can be estimated. For example, the different charm and 
bottom quark mass values which enter into the power law, logarithmic. Coulomb plus lin- 



ear and QCD motivated Buchmiiller-Tye potentials tabulated in ref. [T^] span the ranges 
1.48 GeV < rric < 1.84 GeV and 4.88 GeV < rub < 5.18 GeV. These intervals may be 
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regarded as setting reasonable bounds for the heavy quark mass parameters. The spread 
in values for radial wavefunctions at the origin calculated in these four different potential 
models similarly provides an approximate indication of color-singlet matrix element uncer- 
tainties. Systematic errors which arise from parton distribution functions and higher order 
QCD corrections can also be assessed by performing several fits with different choices of 
distribution functions and renormalization scale. We have not attempted to carry out a 
detailed analysis of the combined impact of all these systematic uncertainties. Our color- 
octet matrix element values therefore represent reasonable estimates rather than precise 
predictions. 

Comparing the numbers in table I with their predecessors in table II of ref. |§] , we see 
that the Q<5['^5'|^^^] matrix element values have all diminished. This is not surprising, for 
some color-octet contributions to quarkonia production are now taken into account by the 
QQ['^Pj^^] and QQ[^Sq^^] channels. We also observe that the NRQCD counting rules are 
more faithfully followed by some matrix elements than others. For instance, the magnitudes 
of {O^^'^i^Si)), {ofi^Si)) and {of {^Po)) /M^ + {of ^ So)) /3 are aU mutually consistent 
with their common scaling rule. On the other hand, {O^"^ {^Si)) is somewhat low while 
{O^^"^ {^Po))/M^ + {of i^So)) /3 is somewhat high. Since ~ 0.23 is not very small, none 
of the charmonia NRQCD order-of-magnitude estimates should be overly interpreted. We 
view the general consistency of the fitted matrix elements with the power counting rules 
as an encouraging indication that the color-octet quarkonia production picture is sound. 

5. Conclusion 

In this article, we have calculated the cross sections for producing colored L = S = 
and L = S = 1 heavy quark-antiquark pairs in hadronic collisions. Intermediate QQ[^Sq^'^] 
and QQ[^Pj^^] states evolve into i/jq mesons at the same order in the NRQCD velocity ex- 
pansion as (5(5['^5'i^''] pairs. Consistency therefore requires that contributions to quarkonia 
production from all three color-octet channels be considered together. We have found that 
the full O(f^) set of color-octet distributions yields good fits to prompt Psi and Upsilon 
data collected at the Tevatron. Numerical values for the long distance matrix elements 
which can be extracted from these data are generally consistent with NRQCD power scaling 
rules. 

Many of the results in this paper can be applied to a range of other interesting prob- 
lems in quarkonium phenomenology. In particular, the NRQCD matrix elements which 
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we have extracted from CDF data are universal and hold for color-octet charmonia and 
bottomonia production at other experimental facilities besides the Tevatron. They can be 
used, for example, to refine the analysis of the J/ij) differential cross section measured at 
the CERN SppS collider which was performed in ref. It would be interesting to see 
whether disparities between gluon fragmentation predictions and UAl data are diminished 
by including QQ\^sf^] and QQ^Pf^ channels The NRQCD matrix elements can 

also be applied to the study of quarkonia production at lepton colliders. Gluon fragmen- 
tation has been shown to represent the largest source of prompt Psi and Upsilon vector 



mesons at LEP ||2^,^. Its incorporation into Z ^ J /ip^ Z ^ and Z ^ T branching 
fractions reduces sizable differences between predictions based upon color-singlet heavy 
quark fragmentation and recent LEP measurements. Color-octet contributions have sim- 



ilarly been found to play an important role in charmonia production at CLEO |25,26 



Finally, the color-octet mechanism may eliminate disagreements between theory and ex- 
periment in fixed target settings. 

We look forward to confronting the color-octet production picture with a variety of 
experimental tests in the near future. @ It will be interesting to see how well this simple 
idea can resolve several problems which currently exist in quarkonium physics. 
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Appendix. Color-octet squared amplitudes 

We list below short distance squared amplitudes for 2 — > 1 and 2 — > 2 scattering 
processes which mediate color-octet quarkonia production. These expressions are averaged 
over initial spins and colors of the two incident partons. The helicity levels of outgoing 
J — 1 and J = 2 pairs are labeled by the subscript h. The total squared amplitudes for 
creating specific quarkonia states are obtained by multiplying these process-independent 
short distance expressions with appropriate long distance NRQCD matrix elements. 

qq ggp^+^Li^^] channel: 

Y,\Am^QWsi'^W = (A.la) 
^|^(,^^QQP^f)])|^=^^ (A.16) 



\h\^i 



gg ^ QQ[2s+ii,(8)] channel: 







])P = 


5(47ra.)2 
192M 


(A.2a) 






])l^ = 


5(47ra,)2 
16M3 


(A.26) 






])P = 





(A.2c) 


h\=i 


ggppf) 







{A.2d) 




ggppf) 


])\' = 





(A.2e) 


E \^^99 ^ 

h\=i 


ggppi^) 


])f = 





(A.2/) 


E \A99^ 

h\=2 


ggppi'^ 


])f = 


(47ra,)2 
12M3 


(A.2^) 



qq QQ[^'+^Lf^]g channel: 



y^lAqq- 












\h\—\ 






E - 






Ei-^(^^- 






E - 

1 1 






Ei^(^^- 

/i=0 




']9)\' 


E 

ihi=i 




']9)? 


E 

\h\=i 




']9)\' 



5(47ras)^ P + 
8(47ra^)3 M'^s 



81M3 (s-M2)4 



81M3 (s - M2)4 

20(47ra,)3 (g - 3i^/2)2(P + ^2) 
81M3 s(s-M2)4 

40(47rQ,)3 s(£2 + ^2) 

81M3 (s - M2)4 

160(47ra^)3 M2t?i 
81M3 (s - M2)4 

8(47ra,)3 g(P ^ ^2^ 
81M3 (s-M2)4 

32(47ras)3 mHu 



[4(P+^2) 



27M3 (s - M2)4 
27M3 s(s - M^y 



(A.3a) 
(A.36) 

£u] (A.3c) 

(A.3(i) 
(A.3e) 

(A.3/) 

(A.3^) 

(A.3/i) 

(A.3i) 
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QQ[^^+^Lf]q channel: 



^ 54M3 [(s-M2)(t-M2)]^ 



h=0 



|h|=i 



X [4{s^ + u^) - su] (A.4c) 



/l=0 



1^1=1 



|ft|=i 



|/i|=2 



27M3 (£-M2)4 
X + ^2 ^ i2M2.«2 ] ( 



9M3 i{i-M^Y 
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gg -> QQ[^''+^L^f]g channel: i 



siu[{s- M^){i- M^){u- M^)Y 



^ 27isi + iu + us)-19M^ 
[{s-M^){i-M^){u-M^)Y 



J: \Agg ^ QWs[%r = -^7-4^ [(^ - + + + 2M^(^)^] 

|7.| = 1 ^ ' 

27{si + iu + us) -19M^ 



[(s-M2)(£-M2)(w-M2)]- 



i7^r3p(8)i.M2_ 5(4™ 



3 



Ei^(^^-wp^on^7)P- 12M3 

X {5254(52 _ ^2)4 ^ ^2^^2(^2 _ ^2)2(3^2 _ 2z^)i2s^ - GsH^ + 3z^) 

+ _ g4^io^2 ^ 2655^5^ - 3825^^6 + 276sH' - 88 fz'^ + 9z''] 

- M^s [54s^° - 3575^^2 + 8445^5^ - 8985^5^ + 4395^5^ - 81z^°] 

+ [I53s^° - 7985^^2 + 14155^^4 - 10415^^^ + 3015^5^ - 18z^°] 

- M^^s [270s^ - 10895^^2 + 13655^54 - QlQfz^ + 87 
+ [324s^ - 9515^^2 + 7695^5'* - 1895^5^ + 9z^] 

- 9M^^s{Qf - e){hs^ - 9ff + 35^) 

+ 2M^^s\brs^ - 5952^2 + 12^4) _ 27M^^s\2f - z^) + QM^'^s'^] 
/[sf{s-M^f{sM'^ + z^f] {k.hd) 



The gg QQ[^ Pj^^] 9 squared amphtudes are expressed in terms of the variables s and 
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h=0 

X sz^ [{s^ - z^f - 2MHz^ - M\f + 2z^) + M^] 



X [{f -z^f- M^s{2f - + mH^] I {{s - M^isM^ + z^] (A.5e) 



■3d(8)i.m2 _ 5(4™ 



\h\=i 



sf 



.6. 



X M2|2(s2 - z^f{s^ - As^z^ + fz^ - 

- MH{2s^ - z^){5f - 17 s^z^ + 95254 _ 56) 
+ M\21s^ - m^z'' + 21sH'' - uH^ + i«) 

+ M^f{lQs^ - 95252 + 2^4) _ M^H\<of - z^) + 
/[(5-M2)4(sM2+52)4] (A.5/) 



3 



/i=0 

+ {S^ - 125652 ^ 4^^4^4 _ 20^2^6 ^ ^8] 

- m65[456- 265^52 _ ^2^4 _ 5^6] 

+ [29s^ - 1145^52 ^ 1085254 _ 1056] 

- M105 [655^ - 1045252 - 335^] + Mi2 [545^ - 205252 + 7i^] 

- Mi45(2352 + 552) + 7Mi652|/ [(5 - m2)6(5M2 + 52)^] (A.5^) 
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\A99^QQ['P^^]9)\ 



X M 




- M6s[42s^ - 7756^2 + 4is4_g4 _ 22 + 17^8] 
+ [53s^ - 8836^2 + 69s^£^ - 68s^z^ + 3z^] 

- M^°s [5456 - 85s^£^ + 60s^£^ - 9^6] + M^^s^ [43s^ - 47s^£^ + 20i^] 

- Ml4s3^22s2 - 9^2) + 5Ml6^4|y _ ^2^6(^^^2 ^ .2^)4j ^^ ^^^ 






306s^^z^ + 4825^5^ - 


2715656 + 77545^ 


- 185251° 




- 5735^5^ + 710^654 ■ 


- 344^456 + 915^5^ 


' - 185^6] 


+ [226s^° - 


- 7415^5^ + 7375654 _ 


310^456 + 77^258 


- 4510] 


-M1°s[264s^ 


- 6865652 + 541^454 • 


- 1775256 + 2558] 




+ M^^[226f - 


- 452^6^2 + 261^45^ - 


555256 + 258] 




-M^^s [I40s6 


- 2015^5^ + 715^5^ - 


656] 





+ MI652 [6O54 - 535252 + 854] _ 2^1853 [852 - 352] + 2M2054} 

/[552(5 -M2)6(5M2+52)4] (A.5«) 




2M3 



X M- 
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Figure Captions 

Fig. 1. Lowest order Feynman graph which mediates q + q Q + Q scattering. 

Fig. 2. Lowest order Feynman graphs which mediate g + g ^ Q + Q scattering. 

Fig. 3. Feynman diagrams which mediate (a) qq i^Qg, (b) gq Vq"? and (c) gg Vgfi' 
scattering through intermediate QQ[^Pf''^] and QQi^S^^^] pairs. 

Fig. 4. Ratio R{p±) of the total QQi^Pf'^] and QQ[^Sq^^] contributions to the tpQ trans- 
verse momentum differential cross section in the limit where the long distance 
NRQCD matrix element {Ot'^i^Po)) equals M^{Ot'^ Sq)). The sohd and dot- 
ted curves illustrate R{p±) for the charmonia and bottomonia sectors respectively. 

Fig. 5. Theoretical transverse momentum differential cross section for prompt tl^' produc- 
tion at the Tevatron in the pseudorapidity interval \r]\ < 0.6 compared against 
preliminary CDF data. The dashed curve depicts the direct color-singlet contribu- 
tion to t/j' production. The dot-dashed curve illustrates the cc[^5'|^^^] cross section, 
and the dotted curve denotes the combined cc[^Pj^^] and cc[^Sq'^] distributions. 
The solid curve equals the sum of the color- singlet and color-octet contributions 
and represents the total theoretical prediction for the differential cross section. 
All curves are multiplied by the muon branching fraction Br('i/'' 

Fig. 6. Theoretical transverse momentum differential cross section for prompt J/i(j pro- 
duction at the Tevatron in the pseudorapidity interval \'r]\ < 0.6 compared against 
preliminary CDF data. The curves in this figure are labeled in the same way 
as those in fig. ^. All curves are multiplied by the muon branching fraction 

Br{J/il; 

Fig. 7. Theoretical transverse momentum differential cross section for J/i/j production 
at the Tevatron in the pseudorapidity interval < 0.6 resulting from radiative 
Xcj decay compared against preliminary CDF data. The dashed curve depicts 
the color-singlet contribution, the dot-dashed curve illustrates the ccf'^S'^^''] cross 
section and the solid curve represents their sum. All curves are multiplied by the 
muon branching fraction Br (J/?/' 

Fig. 8. Theoretical transverse momentum differential cross section for T(IS') production 
at the Tevatron in the rapidity interval \y\ < 0.4 compared against preliminary 
CDF data. The dashed curve depicts the color-singlet contribution which includes 
direct T(IS') production as well as radiative feeddown from X6j(l-P) and XbjC^P) 
states. The dot-dashed curve illustrates the cross section, and the dotted 

curve denotes the combined bb[^Pj^^] and bb[^SQ^^] distributions. The solid curve 
equals the sum of the color-singlet and color-octet contributions and represents 
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the total theoretical prediction for the T(15') differential cross section. All curves 
are multiplied by the muon branching fraction Br(T(lS') 

Fig. 9. Theoretical transverse momentum differential cross section for T(2S) production 
at the Tevatron in the rapidity interval \y\ < 0.4 compared against preliminary 
CDF data. The curves in this figure are labeled the same as those in fig. ^ 
The dashed color-singlet cross section includes T(25') production and radiative 
feeddown from Xbji'^P)- AH curves are multiplied by the muon branching fraction 
Bt{T{2S) ^ ix+ix-). 
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